Stability of the pentaquark in a naive string model 
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The pentaquark is studied in a simple model of confinement where the quarks 
and the antiquark are linked by flux tubes of minimal cumulated length, and the 
Coulomb-like interaction, the spin-dependent terms and the antisymmetrization con- 
straints are neglected.. The ground-state is found to be stable against spontaneous 
dissociation into a meson and a baryon, both in the case of five equal-mass con- 
stituents and for a static quark or antiquark surrounded by four equal masses. 

PACS numbers: 12.39.Mk,12.39.Jh,12.38.Aw 



I. INTRODUCTION 

In the 60s, some indications came up for a possible "Z" baryon resonance with strangeness 
5" = +1, but the data on the kaon-nucleon interaction did not confirm the existence of this 
state. For references, see early issues of the Review of Particle Properties, as cited in the 
last one [1] . In the modern language, such a Z resonance would be a {snnnn) state, where 
n denotes a light quark u or d. 

In the 70s, some excitement arose on baryonium candidates, tentatively interpreted as 
{nn — nn) states with a separation between the {nn) diquark and the {nn) antidiquark, 
and possibly an exotic color charge for these diquarks, see, e.g., [2]. It was then suggested 
that multiquark baryons could exist as well, with a structure {qq — qqq), and again, a 
possible orbital barrier between the clusters and perhaps a color-octet content of each cluster, 
this preventing from immediate rearrangement and subsequent decay into two color-singlet 
hadrons [3, 4]. This work was abandoned when the evidence for baryonium faded away. 

In 1977, JafTe suggested that the dibaryon H{uuddss) could be bound below the lowest 
threshold K{uds) + h.{uds) due to a coherence in the chromomagnetic interaction [5]. In 1987, 
Lipkin [6] and, independently, Gignoux et al. [7], pointed out that the same mechanism 
would bind a heavy pentaquark (the word was invented in these circumstances) such as 
P = (cuuds), (cudds) or {cudss). This pentaquark was searched for in an experiment at 
Fermilab [8], which turned out not conclusive. Further works indicated that the stability of 
H and P hardly survives a more consistent treatment of the short-range correlations which 
enter the chromomagnetic matrix elements and the breaking of the SU(3) flavor symmetry 
in the light quark sector [9-11]. 

More recently, it was shown [12] that in some models of chiral dynamics, new baryons 
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are predicted, in particular an antidecuplet above the usual octet (iV, A, . . .) and decuplet 
(A, . . . , Q^). This triggered a search by Nakano et al. (LEPS collaboration) [13], who found 
evidence for a baryon with strangeness S = +1. Much confusion followed, as stressed in [1], 
where a critical review can be found, with skeptical remarks not so much on the pioneering 
theoretical speculation and experimental search, but on the followers. See, also, [14]. Indeed, 
custom models were quickly designed, where the pentaquark was found with either positive 
and negative parity, made of ad-hoc quark clusters. Several experimentalists discovered 
the potential of their set-up and stored data for looking at exotics and hastily constructed 
mass spectra that could have been investigated much earlier. Eventually, experiments with 
high statistics and good particle identification found no confirmation of the pentaquark 
candidates [1]. There is now reasonable consensus that the light pentaquark does not exist, 
though some puzzling positive indications are still reported [15]. For a recent discussion, 
see, e.g., [16]. 

Nevertheless, the question of multiquarks remains important. On the experimental side, 
several states have been discovered in the hidden-charm sector [1], whose properties suggest 
large {cqcq) components, where q is light or strange. On the theoretical sides, multiquark 
states are now studied with the QCD sum rules [17, 18] and the lattice QCD [19]. In the 
past, the issue of multiquark states was mostly the field of constituent models: in principle, 
the basic ingredients can be tuned by fitting the spectrum of ordinary mesons and baryons, 
and then applied to tentative multiquark configurations. But the main difficulty lies in 
extrapolating the potential from the meson sector to larger systems. For the Coulomb-like 
part, in particular one-gluon exchange, the color additive rule 

i<j 

is probably justified. Here, A^*^^ denotes the 8-vector color generator of the i^^ quark (with 
a suitable change for the antiquark) and the normalization is such that v{r) holds for a 
color-singlet quark-antiquark meson. 

However, there is no reason to use additive rule (1) for the confining part, although it 
has been often adopted as a tentative approximation. For baryons, a y-shape interaction 
has been suggested years ago [20] , and often rediscovered either in models or in attempts to 
solve the QCD in the strong coupling limit [21-26]. This interaction is now confirmed by 
lattice QCD [27] (for further refs., see., e.g., [28]). It reads 

3 

Vy = a min ^ ru , (2) 
1=1 

where each quark is linked to a junction J whose location is optimized, as in the famous 
problem of Fermat and Torricelli. See Fig. 1. Instead, the rule (1) with v{r) — ar would 
give a potential 

Va = 2 (^12 + ^22, + rsi) , (3) 

which is smaller but very close to Vy, so that the phenomenology hardly distinguishes 

between the two potentials Vy and V3 [29]. 

In the case of tetraquarks (two quarks and two antiquarks), however, it was shown that the 
color-additive rule (1) and the generalization of the F-shape potential lead to rather different 
spectra. In the former case, the stability of the tetraquarks {qqQQ) requires a large mass 
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ratio for the quarks and antiquarks. The latter potential, if alone and acting without any 
antisymmetrization constraint (e.g., with quarks and antiquarks of different flavors), gives 
stable tetraquarks [30] . The tetraquark potential was taken as the minimum energy of two 
separate quark-antiquark flux tubes (the so-called "flip-flop" interaction) and a connected 
double-y Steiner tree linking the quarks to the antiquarks, a model of confinement which is 
supported by lattice QCD [31]. 

Our aim is to extend the study of Vijande et al. to pentaquark states. For simplicity, 
we assume that the constituents have the same mass, but remain distinguishable through 
their spin and flavor degrees of freedom, so that the orbital wave function can contain a 
symmetric component. We also consider the case where one of the constituents is inflnitely 
massive. This paper is organized as follows. The model is described in Sec. II. The results 
are shown in Sec. Ill, and some further investigations are suggested in Sec. IV. 



II. THE MODEL 

We focus on the role of conflning forces, and hence disregard the Coulomb-hke contribu- 
tions and spin-dependent forces. For (qq) mesons, the Hamiltonian of the relative motion 
reads 

H2 = — + ar, (4) 
m 

where p is conjugate to the quark-antiquark separation r, r = |r|, and m is the constituent 
mass. For this system and the ones considered below, it is possible to set m = o" = 1 
without loss of generality, since departing from these values results in a simple scale factor 
j7j,-i/3(j2/3 g^Yl the eigenvalues. The ground-state (qq) of (4) can be expressed in terms of 
the Airy function, and its energy is E2 ~ 2.33811. For (Qq) or (qQ) with a static quark or 
antiquark, and a constituent of mass m = 1, the reduced mass is twice larger, and by scaling 
the ground-state energy is E2 ^ 1.8558 

For (qqq) baryons, we consider first the additive model (3), 

H3^pl+pl + ^{ri2 + r23 + rsi) , (5) 

where x — r2 — Vi and y — (27*3 — r2 — ri)/\/3 are Jacobi variables suited for equal 
masses, and and Py their conjugate momenta. We also study the more realistic F-shape 
interaction (2), schematically pictured in Fig. 1, with the Hamiltonian 

Hy^pI+pI + Mri, r2, rg) . (6) 

For (Qqq), the kinetic-energy part is replaced by {pf -|-p|)/2. 

For the five-body problem with equal masses, we start, as in [32], from a symmetrized 
model 

i=l i<j 

where the cumulated strength encountered in a meson and in a baryon within the additive 
model is spread over the 10 interacting pairs. Here, the relative motion is described by four 
Jacobi variables cCj, to be specified shortly, and the p^ are their conjugate momenta. Up to 
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FIG. 1: (Coloer on line) Confinement of mesons and baryons. The minimum over the quark 
permutations gives the flip-flop potential. 

an irrelevant factor, the potential in (7) is identical to the confining interaction adopted in 
Ref. [33]. 

Note that, in the additive model, the threshold can be understood as another five-body 
problem, 

4 

^th = 5Zp- +(l-e)ri2 + 

i=l 

in the limit where e — 0. Then the variational principle applied to Hth with the symmetric 
ground-state of as trial function immediately indicates that the lowest energy of is 
above the threshold. In other words, in the color-additive model (1), the most asymmetric 
distribution of couplings is encountered in a threshold made of two separate color singlets, 
and this asymmetry benefits the threshold and penalizes tentative multiquarks. This is 
confirmed by the results of Hiyama et al. [34, 35]. To build stable multiquarks, one can 
either introduce a competing asymmetry by using different constituent masses, as done for 
the tetraquark (QQqq), or modify the color-additive model [36]. 

For the pentaquark, we consider the natural extension of the minimal-path model already 
used for the tetraquark [30, 37, 38] and supported by lattice studies [31]. It reads 

4 

Hp = Yl P' + ^P^ ^P = min(^ff' ^st) , (9) 

where 

Vff = min [rii + Vyirj, r^, Vg)] , (10) 

i 

with {i,j,k,i} being any permutation of {2,3,4,5}, is the so-called flip-flop potential. It 
corresponds to the most economical configuration in Fig. 1. The second term, Vst, cor- 
responds to a connected Steiner tree, as pictured in Fig. 2. This Steiner tree generalizes 
the Fermat-Torricelli problem for more than three points. The quarks {i,j) are linked to 
the central Steiner point S as three quarks in an ordinary baryon, through an intermediate 
Steiner point Sij. Similarly, {S,k,£) form a baryon-like structure with Steiner point Ski- 
Then, the antiquark and the two intermediate Steiner points Sij and Ski form an antibaryon- 
like configuration. Of course, the potential Vst is optimized by varying the permutation 
{i,j,k,£} of the quarks. See, e.g., [38, 39] for references on the Steiner problem and its 
application to the multiquark potential. 

Note that the long-range part of the pentaquark potential, as estimated in lattice QCD, 
has been found fully compatible with this multi-F term [40]. 
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FIG. 2: (Color on line) Example of connected Steiner tree (planar for simplicity) linking the 
antiquark 1 to the quarks {i,j,l,£}. The Steiner points are S, Sij and s^e- 



An obvious consequence of the minimization of (9) is that 



Vp < ri2 + VV(r3,r4,r5 



(11) 



indicating that the pentaquark potential is smaller than the cumulated confinement energy 
of the threshold, i.e., that the effective potential between the meson (1,2) and the baryon 
(3,4,5) is attractive. However, in three space dimensions, an attractive potential does not 
automatically lead to binding, if this potential is short-ranged. Hence one should solve the 
five-body problem to determine whether this model supports or not a stable pentaquark. 

In the case of the tetraquark, it was noticed [30] that the connected Steiner tree config- 
uration, though the most interesting, plays a minor role, and that the binding is obtained 
from the fiip-fiop term alone. Similarly, a survey with randomly generated coordinates for 
the antiquark and the four quarks shows that, in Eq. (9), the Steiner tree gives the minimum 
in less than 3% of the cases. Hence, we shall neglect this term and compute an upper bound 
with the fiip-fiop interaction alone. This is opposite to the choice made in [32], where only 
the connected-Steiner-tree term is adopted and the fiip-fiop one omitted. 



or, in short, |\E'2) = Ylilil'^i)^ ^^'^ t^e energy E2 — 2.33811 can be reproduced with 
just a few terms. The relevant matrix elements {a'\a), {a'\p'^\a) and (Q;'|r|Q;) are known 
analytically, and are basic ingredients for more complicated systems. 
For the ground-state of if 3, a generalization reads 



III. RESULTS 



The ground-state meson can be described by a simple expansion 




(12) 




(13) 



and, again, the matrix elements are known analytically. For a given choice of range parame- 
ters, the weights 7^ are given by a generalized eigenvalue equation. To avoid ambiguities and 
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simplify the minimization, one can restrict the Gaussians to scalar (a = 6, c = 0) or diago- 
nal {a ^ b, c = 0) matrices, and those given by permutation of the quarks; furthermore the 
parameters Oj and 6j can be taken from a single set {a, a + 6,a + 25, . . .}, with minimization 
only over the two extreme values [41]. One reaches ~ 3.863. If only scalar matrices are 
allowed, the expansion (13) converges towards the best function of the hypcrradius given by 
p2 — 3.2 _|_ ^2 ^Yiis approximation, u{p) — p^/^ ^ is given by the radial equation 



-u"{p) + ^«(p) + Voopuip) = EfK{p) , (14) 



with suitable boundary conditions, which leads to the upper bound Ef^ ~ 3.865. The 
hyperscalar projection, including the strength factor 1/2 and the number of pairs, is 

T/qo = ^ ^ 1.019 . (15) 

OTT 

A similar strategy can be used for the F-shape potential, except that the matrix elements 
have to be calculated numerically. The hyperscalar coefficient becomes Vqq = 1.115. By 
scahng from (15), this corresponds to an energy E^^ ~ 4.105. The Gaussian expansion, if 
not restricted to scalar matrices, gives a better energy Ey — 4.095. This means that, in our 
simple string model, the threshold for the stability of light pentaquark states is 

Et^^iqqqqq) ^ 6.433 . (16) 

If a static quark Q or antiquark Q is introduced, for (Qqqqq), the threshold consists of 
{Qq) + {qqq)j ^^^^ i^i {qqqqQ), it is the lowest of {qQ) + {qqq) and {qq) + {Qqq), which turns 
out to be the latter by a small margin, and the thresholds are 

Et^Qqqqq) ^ 5.950 , E,^{qqqqQ) ~ 5.944 . (17) 

We now turn to the five-body problem, first for (qqqqq). A possible choice of Jacobi variables, 
besides the center-of-mass, is 

4ri - (r2 + r3 + r4 + rs) ^2 - r4 

VlO 2 .^gx 

_r2 + r3-r4 + r5 _r2 + r3 + r4-r5 

X3 - - , X4 - - , 

which is convenient to express the permutations of the quarks. With a Gaussian expansion, 
one obtains for the symmetric toy model an energy Er, ~ 6.850. In the hyperscalar 
approximation, the radial equation is similar to (14) with now a centrifugal coefficient 99/4 
instead of 15/4, and a strength Vqo — 2560/(6937r) ~ 1.176, and a very similar energy is 
found. This model has thus energies well above the threshold E2 + -Es that includes 
pairwise forces for the baryon, and even above the threshold E^]^ = E2 + Ey when the 
baryon is bound by the F-shape interaction. 

If now one switches to the flip-flop interaction, the coefficient of the hyperscalar potential 
Vqo P is found to be Vqo = 1.031, by numerical integration of Vg over the hyperscalar variables, 
for given p. By scaling, this gives an energy 

E^g\qqqqq) 6.276 , (19) 
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in the hyperscalar approximation, which is clearly below the dissociation threshold (16). 
This is confirmed by the Gaussian expansion. As these two methods are variational, the 
flip-flop model gives a bound state, and this is a fortiori the case if the connected Steiner 
trees are also included in evaluating the potential (9). 

If the calculation is repeated in the case of four unit masses and an infinitely massive 
antiquark or quark, the kinetic energy in (9) is replaced by Ylpff^^ where is conjugate 
to the position Vi of a finite-mass constituent. The pairwise model (7) gives an energy of 
about 6.305 for both {Qqqqq) and {qqqqQ) in the hyperscalar approximation, i.e., well above 
the threshold, and also with the Gaussian expansion, no state is found below the threshold. 
On the other hand, the flip-flop potential, treated in the hyperscalar approximation, gives 

eP (Qqqqq) ~ 5.836 , (qqqqQ) ~ 5.667 , (20) 

which is sufficient to demonstrate the of stability of (Qqqqq) and (qqqqQ) with respect to 
their respective thresholds (17). 

IV. CONCLUSIONS AND OUTLOOK 

A simple string model of linear confinement gives a pentaquark which is stable against 
spontaneous dissociation into a meson and an isolated baryon. This is at variance with most 
of the earlier constituent-model calculations. There are, however, severe limitations in our 
approach: 

• The non-relativistic kinematics for the quarks and the Born-Oppenheimer treatment 
of the gluon field, which supposedly readjusts itself immediately when the constituent 
move, call for an application to heavy quarks. However, the short-range central cor- 
rections should be incorporated for heavy quarks. For the Coulomb-like interaction, 
the color- additive rule (1) probably holds. If alone, a Coulomb interaction with color 
factors will not bind. Our pure confining model binds. What occurs for a superpo- 
sition would be a matter of detailed phenomenology beyond the scope of the present 
note. 

• The quark wave function is assumed to be compatible with an overall s-wave, and is 
symmetric under permutations. Thus, it should be associated with enough spin and 

fiavor degrees of freedom in the quark sector. A proper treatment of the Fermi statistics 
of quarks is rather delicate in the fiip-flop or Stcincr-trcc model, as different flux tube 
topologies correspond to different color couplings, and thus different constraints for 
the spin, flavor and space parts of the wave-function. 

• The calculation has been restricted to flve equal and flnite masses, or one inflnitely 
massive constituent surrounded by four equal masses. The property of stability was 
found to survive for both (Qqqqq) and (qqqqQ). The investigation should be extended 
with better variational wave functions, and a larger variety of mass distributions for 
the constituents. Clearly the conflguration (cuuds) and its analogs by permuting the 
light quarks or by replacing c by 6 would deserve a reflned treatment including the 
spin-spin forces, which are favorable [6, 7]. 

• The pentaquark states with two or more heavy constituents would deserve a speciflc 
studies. For (QQq) baryons, the dynamics can be studied in the Born-Oppenheimer 
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approximation, where the two heavy quarks experience an effective interaction re- 
sulting from their direct interaction modified by the hght quark [42-44]. Similarly, 
the (QQqq) tetraquarks and the (QQqqq) or (qqqQQ) pentaquarks could be studied 
through in the adiabatic limit. 

It remains that our conclusion, based on the flip-flop dynamics, is drastically different 
from the one obtained from the connected Steiner tree alone [32] . Studies within the lattice 
QCD [31] or, more recently, the AdS/QCD [45], have analyzed the interplay between flip- 
flop and connected multi-F conflgurations for tetraquarks. It would be desirable to have 
more information about the analogue for pentaquarks. 
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